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Abstract 

In this article the statistical properties of symmetrical random matrices whose elements are 
drawn from a g-parametrized non-extensive statistics power-law distribution are investigated. In 
the limit as g — > 1 the well known Gaussian orthogonal ensemble (GOE) results are recovered. The 
relevant level spacing distribution is derived and one obtains a suitably generalized non-extensive 
Wigner distribution which depends on the value of the tunable non-extensivity parameter q. This 
non-extensive Wigner distribution can be seen to be a one-parameter level-spacing distribution 
that allows one to interpolate between chaotic and nearly integrable regimes. 
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The Hamiltonian M x M matrix that is examined is symmetric and real. In order to 
simphfy the derivation M is restricted to M = 2. The derivation of the 2x2 random 
matrix statistics follows the pedagogical approach of Brody a/ 0|. A discussion of 
random matrix theory as applied to quantum chaos and integrability is found in [0, The 
Hamiltonian matrix is 



Hi2 H22 

The first task in random matrix theory is to obtain the probability of occurrence of the 
individual elements in the ensemble. This is obtained in this derivation from the maximum 
entropy approach p. Also, in order to guide the subsequent non-extensive statistics deriva- 
tion, the Gaussian distributed random matrix results are derived first and then generalized 
to the non-extensive case. 

The question of the probability of occurrence of the individual elements of the Hamilto- 
nian matrix in the maximum entropy approach is one of obtaining the central moments, i.e. 
the means, variances of the matrix elements. In order to obtain the most likely (least biased) 
probability density for the occurrence of the matrix elements the entropy of the system is 

N 

maximized (S) given the set of constraints (the observables) J2 i^i) 

i=l 



(S) = - J P{Hu,H22,H,2) lnP{Hn,H22,H^2)dHudH22dH,2, 

(-f^ii + H22 + 2^/12) ) ^ / (-^11 + -^22 + 2-^^12) PiHii, H22, Hi2)dHiidH22dHi2 

(2) 

The maximization of the entropy is 

6 {S) - 6[(3{{Hl, + Hi, + 2Hl,))] ^ 0, (3) 

and here is a Lagrange multiplier which is related to the variance by /? = The 
maximization yields the least biased probability density [A is the normalization) 

P{H) = P{Hn,H22,Hu) 

= A e"^(^ii+^22+^^"). (4) 



The probability density is seen to be a Gaussian and satisfies the normahzation condition 

J P{H)dHndH22dHi2 = 1. (5) 

Furthermore, this form of distribution has the following properties that are of importance 
in the random matrix theory. 

(i) The probability density of the statistically independent random matrix elements is 
factorizable. This is written as 

P(Hn, H22, Hu) = P(Hu)P(H22)P(Hu), (6) 

and can be seen by inspection to be true. This property also indicates that there are no 
correlations in the probabilty density as the individual elements are perforce statistically 
independent. 

(ii) There is a related issue of extensivity (additivity). The entropy measure for the system 
is formed from the addition of the entropies of the parts such that (omitting constants) 

S = lnPiHn,H22,H,2)^ln[PiH,,)P{H22)PiH,2)] 

= lnP{Hn) + lnP{H22) + lnP{Hu). (7) 

(iii) The probability density is invariant under orthogonal transformations of the form 



H' = O'^HO, 

cos {9) sin{9) 
—sin{6) cos{6) 



O = 



(8) 



This transformation yields the following transformed variables in the infinitesimal limit of 

i^n = Hii — 26Hi2 
= H22 + '29Hi2 

= H^2 + 9{H^,-H22). (9) 
A transformation of the probability density is written as 

P{H)dH = P{H')dH', (10) 
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and a calculation of the determinant of the Jacobian of the orthogonal transformation = 
det{J) = 1 yields 

P{H)=P{H'). (11) 

This probability density is invariant under orthogonal transformations, and the level statis- 
tics obtained is known as The Gaussian orthogonal ensemble (GOE). Unitary and simplectic 
transformations and symmetry considerations yield the Gaussian unitary ensemble (GUE) 
and the Gaussian simplectic ensemble (GSE) in a related fashion, however these cases of 
symmetries and degeneracies will not be discussed further in this article. 

The extensive Gaussian random matrix theory can be generalized by examining the non- 
extensive, or g-parametrized entropy |^, § (equivalently, an incomplete information mea- 
sure). For systems with statistically dependent (say, Hamiltonian matrix elements) variables 
the joint probability decomposition is 

P{H,,H,) = P{H,\H,)P{H,), (12) 

which gives the pseudo-additive entropy 

S,{H,,H,) = Sg{H,) + Sg{H, I H,) + (1 - q)S,{H,)Sg{H, \ H,), 

pl-q _ I 

S, = -ln,P = . (13) 

It is known that the Tsallis entropy satisfies this condition, and the resulting probability 
will be of the power-law form. The g-logarithm is IrigX = — ^-^^yy-. In the limit as g — 1 
the usual form of the natural logarithm and thus the extensive statistics and its exponential 
(Gaussian) distributions is recovered. 

The entropy to be maximized given the constraints is then 



J [Hl, + Hl^ + 2Hl^) P'^{Hn,H22,H,2)dHudH22dH^2 

(14) 



and which is subject to the extra normalization condition 



P{Hu, H22, H^2)dHndH22dH,2 = 1. (15) 
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The maximization is then similar to the Gaussian case, and the variation given the 
constraints is 

6 (S)^ - 6[f3{{Hl + Hi, + 2Hl,)) ] ^ 0, (16) 



which upon solving yields the least biased probability density 
P(iJn, Hu) =A (l+l3{q- l){Hl^ + iJ| 
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(17) 



This is then the Tsallis power-law form of the probability density for the 2x2 random 
matrix elements. Having obtained the statistics of the random matrix, the explicit form of 
the non-extensive Wigner distribution is derived. This non-extensive Wigner distribution 
will be the level spacing statistics for a system with statistically dependent random matrix 
elements. 

The eigenvalues of the 2x2 Hamiltonian matrix, Eq.(|l]), are given by 
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2 ' ' 2 

which can be written in terms of a diagonal matrix 



H. 



22 



^-"12 



(18) 



D 



E+ 
E_ 



(19) 



(20) 



This matrix is related to the Hamiltonian matrix by another orthogonal transformation 
H = non^ where is 

cos(0) —sin{(j)) 
sin{(f)) cos{(f)) 
This transformation then relates the variables as 

Hii = E+cos^{(t)) + E_sin^{(j)) 
H22 = E^sin 

H12 = {E+ - E^)cos{(f))sin{(j)). 



(21) 



Next the probability densities are transformed directly as 

P{Hu, H22, Hu)dHndH22dHu = P{E+, E_, (f))dE+dE_d(f), 

and calculating the determinant of the Jacobian of the transformation gives 

d{Hii, H22, H12) 



det{J) = det- 



d{E+,E_,, 



E^-E_ . 



(22) 



(23) 



Rewriting the argument of the probabihty density in terms of the new variables 

^11 + ^22 + 2i^i'2 (24) 

shows that the probabihty density is clearly independent of . The transformed probability 
can then be written as {P{E+, E^, cj)) = P{H)det{J)) 



P{E+, E_) =A {E+- E_) (l + (5{q - l){El + El 
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(25) 



In order to obtain the non-extensive Wigner distribntion, the variables arc recast in terms 
of the level spacing and an auxiliary variable that plays the role of a 'center of mass' energy 
coordinate 



s = E+-E_ 
E+ + E. 



z — 



(26) 



Substitution of these variables into the probability density results in 



P{s,z)^As l + -!)(_ + 2^2) 
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(27) 



and integration over the auxiliary variable z yields the level spacing distribution 
P{s) = A 



V 2^ 


1 i' 

q-l 2 




2^/?(g-l)r 


1 



S l+/3(g-l)(-) 



, -1 I 1 
-,2 \ -^+2 



(28) 



Next the level spacing distribution P{s) is normalized to a mean level spacing of one (s) = 1 
and one obtains the normalized Wigner distribution 



P{s)^^(5{q-l) y-^' l{ s (l + /?(g-l).V2)'-"% (29) 
/ t(— 2) \ ^ 

where is given by /3 (g) = (^ r(^'-f ) j ' 

The behavior of the non-extensive Wigner distribution is obtained by plotting its values 
for a range of the level spacing s, given the 'inverse variance' 13 and the nonextensivity 
parameter q. In Fig.l. the q dependence of (3 is plotted for values of q between 1 < q < 1.4. 
In Fig. 2. the Poisson (sohd), extensive (long-short dashing) and non-extensive Wigner (short 
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FIG. 1: /3 Vs. q 

dashing) distributions are plotted for a low non-extensivity parameter q value of g = 1.01. 
The non-extensive distribution is nearly superimposed on the extensive Wigner distribution 
as is expected for q— > 1. In Fig.3. the Poisson (solid), extensive (long-short dashing) and 

non-extensive Wigner (short dashing) arc plotted for a high value of the non-cxtcnsivity 
parameter q = 1.38. Here the distribution is greatly shifted and approaches the Poisson 
level statistics distribution. 



In this letter the Gaussian orthogonal ensemble (GOE) results for a 2 x 2 random Hamilto- 
nian matrix is generalized to the case of the non-extensive statistics and the resultant power- 
law distributions. A derivation of the subsequent level spacing statistical distribution, the 
non-extensive Wigner distribution, is given. This derivation is obtained by maximizing the 
Tsallis non-extensive entropy for the 2x2 symmetrical random matrix elements. This can 
be straight-forwardly generalized to M x M matrices. The resultant non-extensive Wigner 
level-spacing distribution is gr-parametrized and allows for a smooth interpolation between 
the extensive Wigner distribution and the regime where the level statistics are given by a 
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FIG. 2: P(s)^ Vs. s, q = 1.01. The non-extensive and extensive Wigner distributions axe nearly 
super-imposed. The Poisson distribution is plotted using a solid line. 

Poisson distribution. In future work it will be interesting to apply these results to Hamil- 
tonians of mixed systems between regular and chaotic regimes where deviations from the 
Wigner statistics become pronounced. 
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